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ABSTRACT 


The mechanism of magnetic pumping consists of two processes, 
the adiabatic motion of charged particles in a time-varying magnetic 
field and their pitch-angle diffusion. The result is a systematic 
increase in the energy of charged particles trapped in mirror (and 
particularly, magneto spheric) magnetic fields. A numerical model of 
the mechanism is constructed, compared with analytic theory where 
possible, and, through elementary exercises, is used to predict the 
consequences of the process for cases that are not tractable by 
analytical means. For energy dependent pitch- angle diffusion rates 
characteristic "two-temperature" distributions are produced. 

The model is applied to the outer Jovian magnetosphere for. 
two purposes; to find magnetospheric regions in which the mechanism 
may energize trapped particles, and to generate distribution func- 
tions involving pitch- angle diffusion caused by wave-particle inter- 
actions. We find that beyond 20 Rj in the outer magnetosphere 
particles may be magnetically pumped to energies of the order of 
1-2 MeV and we predict two- temperature distribution functions with 
"break points" at 1 - 4 KeV for electrons and 8-35 KeV for ions. 
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INTRODUCTION 


Unexpectedly large fluxes of energetic particles were de- 
tected in the outer Jovian magnetosphere by the Pioneer and Voyager 
spacecraft measurements. With energies too high to be explainable 
by solar-wind injection, the source of these particles is assumed to 
be within the magnetosphere. One possibility for the energization, 
proposed by Goertz [1978], is the magnetic -pumping process of Alfven 
[1954] . In Alfvln's model [Alfven and Falthammar, 1965 ] cyclic 
variations of the magnetic field strength are combined with iso- 
tropization of the particle distribution functions when the field 
strength reaches a maximum and/or minimum. As the field strength is 
increased, the mean perpendicular (to 5) momentum of the distribu- 
tion increases, conserving the first adiabatic invariant of each 
particle, while the mean parallel momentum remains constant. Iso- 
tropization of the distribution when the field is at a maximum turns 
some of the high perpendicular momenta into high parallel momenta, 
which do not decrease when the magnetic field is reduced. Repeating 
this cycle leads to an increase in the mean momentum of the distri- 
bution with time. 

Goertz ' s model is based on the idea that low-energy magneto- 
spheric particles, conserving their total (kinetic plus potential) 




energies, will 2 x ^ drift almost precisely on electric equipoten- 
tials. If the magnetic field is not rotationally symmetric, the 
drift will carry the particles through regions of differing magnetic 
field strength, hence the magnetic cycle. Pitch-angle scattering of 
the particles is assumed to he effected by the presence of plasma 
waves, for example, those observed by the Voyager plasma wave 
experiment [Gurnett et al., 1979]. One shortcoming of the Goertz 
model is the requirement that pitch- angle scattering occur only at 
local noon and midnight. Energization through the magnetic pumping 
process will also result from various amounts of pitch-angle scat- 
tering occurring, more realistically, throughout the magnetosphere. 
Furthermore, Goertz did not evaluate steady- state distribution func- 
tions. Both of these shortcomings will be corrected in this paper 
although others still remain. 

We begin in Section 1 by studying the foundations necessary 
to model the process of magnetic pumping. We use the conservation 
of two adiabatic invariants and an assumed shape for the magnetic 
field to follow the adiabatic motion of particles. In addition, we 
solve, at each time step, a diffusion equation describing the pitch- 
angle scattering of particles. 

In Section 2 we combine the two devices, adiabatic motion and 
pitch-angle scattering, to simulate the magnetic -pumping process. 

We first discuss the theory of the extreme, Alfvenic case, and then 
complete several exercises for the more realistic case with 
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pitch-angle diffusion occurring simultaneously with the changing of 
the mirror magnetic fields. We will identify requirements necessary 
to make the particle energization most efficient, look at the time 
hehayior of the mean momentum of a group of particles, and examine 
several distribution functions generated by the mechanism. 

In the application to Jupiter, Section 3, we concentrate on 
predicting the magnetospheric regions in which magnetic pumping may 
be expected to energize charged particles before they are lost, and 
on examining distribution functions produced by diffusion coeffi- 
cients modeled upon resonances with plasma waves. These results 
pertain to observations made by s;*.- see craft in the Jovian magneto- 
sphere of thus far unexplained energetic particle distributions. 
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1. ADIABATIC fOlCN IN MIRROR GEOMETRY AND 
DIFFUSION IN PITCH-ANGLE SPACE 


Using adiabatic theory, our goal is to follow the evolution 
in phase space of a distribution function of particles trapped in a 
magnetic flux tube. 

Changes in the field parameters of a flux tube as it convects 
through a nonsymmetric magnetosphere, or, in the case of mirror 
machines, as the laboratory magnetic field is altered, will produce 
distortions of the distribution. We follow separately each element 
of n(P,o:), the number of particles in a flux tube with momentum P 
and pitch angle a, calculate the new phase space coordinates for 
each element, and reassemble the distribution. The method used will 
not require the field to be known for all values of time. Given a 
distribution in a flux tube with known field parameters, a descrip- 
tion of the new field topology will be sufficient to calculate the 
new distribution. 

The coordinate transformations are obtained, from the two 
adiabatic constants of the motion for mirror geometry. An idealized 
magnetic field topology described below, which although not neces- 
sary, will simplify our calculations. The first adiabatic invariant 
(magnetic moment), may be written 
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(1) 


Using the particle's pitch angle, a (a = 0° for field-aligned 
velocities), we esqcress the perpendicular momentum of a particle as 
= P sin a, yielding 

sin^ a = jp- , (2) 

m 


where is the mirror field strength. With this equation and the 
knowledge of a particle ' s pitch angle at a point of known magnetic 
field strength, the mirror field strength (independent of particle 
momenta) can be determined. 

The second adiabatic invariant (longitudinal invariant), J, 
may be written 

dz , (3) 

o" 


where z is the distance along a field line from the equatorial (B = 
minimum) region. The limits of the integral are thus z = 0 (equa- 
tor) to z = z (mirror point). Writing P as P = P'/ 1 - sin^ a and 
' m ' II ! 

2 

using Equation (2) to eliminate sin a, this becomes 
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z 

J * P J m Jl - B/B dz 
o li 


( 4 ) 


where B = B(z). Dividing J by we obtain a hybrid invariant, 

K = v/b^Tb dz . (5) 

We now assume a model for the magnetic field strength, B(z), 
along a particular flux tube, taking it to be of the form [Goertz, 
1978] 


B = B + \ a z 2 (6) 

o 2 • ' ' 

with the restriction that there exists some z = z such that 

c 

B(z c ) = B q = constant (7) 

(see Figure l) . Although this model cannot be regarded as realistic 
for high latitudes, its analytical convenience is an obvious advan- 
tage. We do not expect the qualitative nature of our results to 
change if we would add higher order terms in Equation (6). We have 
tested the effect of varying the value of z c and B(z c ) and have 
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found no change in energization rates. We also used a square root 
dependence of B(z) in z and have found little difference in the 
energization rates. 

Inserting the magnetic field model (6) into Equation (5) for 
the hybrid invariant K, and using Equation (6) to obtain the upper 
limit of integration, we find 


K •= 


TT 


(B - B ) 

2^2 Jg m ° 


(8) 


We are now ready to follow the adiabatic evolution of a distribution 
function through momentum-pitch- angle space. If the equatorial mag- 
netic field strength is changed from its initial value B^ E to a new 

NEW , . . . . 

value B , the invariance of B(z ) in the field model (6) requires 
o c 


2f OLD _ NEW, 

NEW OLD ''o o 1 

a = a + g . (9) 

z 

c 

Thus, for our field model, with any value of a known, the magnetic 
field for all z is entirely specified by its equatorial value, B q . 

When B q changes, an element in phase space with equatorial 
coordinates P^ 0 , will move to the equatorial coordinates P^ E ^, 

ofW 

Using the conservation of the hybrid invariant, Equation (8), 
the particle's new magnetic mirror field strength is found to be 
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-NEW _NEW /^OLD -OLDv /a*’ 
m o x m ° V 


where (Equation (2)) 


. 2 OLD 
sm o: 

o 


Applying Equation (2’’ in the equatorial region of the new magnetic 
field with Equation (10) yields an equation for the new equatorial 
pitch angle in terms of known quantities, 


(^NEW _ 
o 


arcsm 


-P.NEW / 

D / 

m ' 


The invariance of the magnetic moment, Equation (l), applied 
at the new and old mirror points yields 


„NEW -OLD / m 

= Vb old 

m 


Thus, our .coordinate transformation is obtained as follows. 

/ \ NEW 

Set the new value for B q and use Equation (9) to obtain a . Then, 
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for each point (P°^ D ,a 0Tj ^), evaluate Equations (IQ), (12), and (13) 
in that order. 

The numerical problems of representing the distribution func- 
tion on a finite phase-space grid is discussed in Appendix I. 

The bounce-averaged, pitch- angle diffusion equation for the 
equatorial particle distribution function f (a o ) [Lyons, Thorne, and 
Kennel, 1972] is 


at " T(a )sin a cos a aa 
o o o o 

( 1*0 

where T(a ) is the angular dependence of the bounce period and jO(Q’ q ) 
is a "bounce- averaged" diffusion coefficient. This equation has 
singularities at a Q = 0 and a^ = tt/2 which disappear when we trans- 
form it into a diffusion equation for the total number of particles 
in a flux tube, n(a ), expressed in terms of equatorial parameters. 
If f (a ) is the equatorial distribution function, then 

f (a )2tt sin a da = number of particles in the equa- 

torial region with equatorial 
pitch angles from a Q to 

a + da 

o o 


T(a )sin a cos a fi(a ) -r~ 
o' o o o' da 


( 15 ) 
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The extent of the region a particle of given pitch angle 0C Q and 

momentum p traverses is twice the distance z m along the flux tube 

from the equator (minimum magnetic field) to the mirror point. This 

may be expressed in terms of v n , the component of a particle's mean 

‘ li 

velocity parallel to the magnetic field (parallel to the flux tube), 
as 



t B 



(16) 


where is the bounce period. With the magnetic field model used 

(Equation (6)), we may show that v = tt/4 v ccs a , enabling us to 

“ll ° ° 

write 


2z = 5 Tt, v cos a 
m o B o o 


(17) 


The total number of particles with equatorial pitch angles between 

a and a + da , n(a )da , will be the product of the number of 
o o o’ o' o’ * 

particles in the equatorial region with equatorial pitch angles 
ranging from a to a^ + da^ and the extent of the region they 
occupy. Thus, from relations (15) and (17), we have 


n(a )da = 2rr v f (a )t t . sin a cos a da 
' o' o 8 o o s o' B o 00 


(18) 
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Using this relation to replace f ( 0 ! o ) in Equation (l 4 ) yields 


M%) 

st 


'B 


T(a ) *a 
s o' o 


Ka^sin a 0 cos a o £ j|- 


n(a o ) 


sin Ci cos 


(19) 


2 

where a factor of 2tt v q /8 has been divideu out. We may factor the 
bounce period t b into T(v o )T(a o ), where T(v q ) and T(a Q ) are the 
velocity and equatorial pitch-angle dependences. On the right-hand 
side of our equation, the two factors of T(v ) will cancel, leaving 


T(a ) sin a cos a &{a ) — — 
' o' ' o o s o' Ba 

o 

n(a o ) ", ■ 

* ■ T(a )sin a cos a j 
o' o o J 



( 20 ) 


(This equation may also be derived from first principles by bounce- 
averaging a Fokker-Planck equation for the diffusion in pitch- angle 
space [Borovsky, 1980, Chap. III].) Using the product rule on the 
inner differentiation of (20) and defining 
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H(« o ) = fi(a 0 )T(a )sln « o cos « o ^ ) sin a cos « 

O O 0 0 


‘(“o' Sa"{ los 

o 


T(a )sin a cos a 
' o o o 


) ' 


( 21 ) 


simplifies our diffusion equation to 


a"( a o> a 

at " da 

O 


dn(a )-, 


p 

*< a o> » 1 

L 0 - 

SOU 

0 

H(a o )n(a Q ) 

■ «• 


(22) 


This diffusion equation describes the pitch-angle scattering 
in terms of the function n(a ), the number of particles in a flux 
tube with equatorial pitch angle a . The numerical treatment of 
Equation (22) is deferred to Appendix I. 
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2. THEORY AND RESULTS FOR MIRROR GEOMETRY 


As a first example of the mirror- geometry energization mech- 
anism we consider Alfven's model of "magnetic pumping" [Alfven and 
^althammar, 1963, Sect. 2.7.4]. We follow, in time, a distribution 
of trapped particles as the equatorial magnetic field is varied from 

an initial, minimum value of B , to a maximum value B , and 

°MIN °MAX 

returned back to the value B . The model requires that the dis- 

°MIN 

tribution be diffused to isotropy in pitch-angle space when the 
field is at its maximum and minimum values. The significant param- 
eter is the ratio k, defined as 


°MAX 

1 

°MIN 


(23) 


The ratio of the mean momentum of a distribution after one cycle 

(B,™ -♦ B lirAV -* B.—,..) to the mean initial momentum is 
v MIN MAX MIN' 


JL 

P. 


1 

3 


5 + 2k + 


k 


(24) 
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Note that, in this model, a value of k = 15.44 will double the 
momentum (quadruple the energy, nonrelativistically) of the trapped 
distribution of particles each cycle. In Figure 2 we compare, as a 
function of the magnetic field ratio k, results of computer simula- 
tions of the model to this theoretical prediction for the momenta 

ratio. The simulations used a 46 x 350 (a Q X P) size grid to repre- 

2 

sent the particle distribution, and used the B(z) = B q + l/2 a z 
magnetic field model in following the adiabatic motion (as B q and a 
varied) of the distribution. Attributable to finite grid effects is 
the fact that the simulations resulted in slightly lower rates of 
momentum gain than analytical theory predicts. This is due to the 
boundary conditions requiring the diffusion coefficient and its 
spatial derivative to be zero (see Appendix I) at CC o = 0, n/2. 

Thus, the distribution on the end points of our grid does not dif- 
fuse with the rest of the distribution. Figure 3 displays the 
evolution in time of the function n(P,a Q ) during one cycle of the 
Alfvenic model. 

We now look inco a more applicable, but more complicated, 
model similar to the one proposed by Murty and Varma [1958]. In the 
Alfvenic model, the diffusion is restricted to times when the mag- 
netic field is at its maximum and minimum values. At those times, 
complete isotropization is allowed to occur. In this more realistic 
model, we require the diffusion to occur continuously, and at a con- 
stant rate, as the equatorial field strength changes (sinusoidally) 


IT 


' 

i 

i 

; with time. There arc two limiting values of the diffusion rate 

i (amount of diffusion /cycle) for which there will be no momentum 

gain. The first is the limit of the diffusion rate going to zero. 

t 

■ Throughout the cycle the distribution will change adiabatienlly (no 

| 

scattering to violate the first or second adiabatic invariant), 
returning to its original form at the end of a cycle. The second 

j case is the limit of infinite diffusion, forcing the distribution to 

remain isotropic at all times as the magnetic field varies. Using 
simple arguments we show in Appendix II that, in this limit, the 
energization approaches zero. 

The results of computer simulations of this model are dis- 
played in Figure 4. Note that there is an optimum diffusion rate 
producing a maximum in the momentum gain per orbit for a distribu- 
tion of particles. Murty and Varma describe this maximum as a 
resonance between the cycle of field variation and the diffusion 
rate. In Figure 4 we have used the diffusion time = 0.0738/$ 
where the constant 0.0738 is valid for the numerical parameters of 
our model. This time can be interpreted as the time it would take 
the lowest eigenfunction of the pitch-angle diffusion operator to 
decay to l/e of its initial value for a constant magnetic field. In 
the case of k = 3 (P = 1/2), we find the ratio of the momentum gain 
per orbit (P - F q )/P q of this model to the momentum gain per orbit 
of the Alfven model to be .520 (Murty and Varma predict a ratio 
of .408). 
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The magnetic ratio k (= B , / B ) dependence of this model 

°MAX °MIN 

is displayed in Figure 2. The diffusion coefficient was chosen so 
that = ’’’CYCLE' S :i - v: ^ n S the maximum obtainable momentum gain per 

cycle . 

Of interest in this model is the mean momentum of a particle 
distribution as a function of time. If, as in Equation (24), we 
find that after one orbit the mean momentum of the distribution 
changes according to the relation 

P 

— = constant , (25) 

o 

then the average value of P as a function of time will be given by 
[see Alfven and Falthammer, 1963] 

P = P e t / T (26) 

o ' ' 

where t = t(B /B ) is time- and momentum- independent. In 
°MAX °MIN 

Figure 5 the results of numerical simulations, for two values of 

k = B /B , are plotted. The envelope of the k = 3 curve fits 
°MAX °MIN 

the form of Equation (26) quite well, with P q = 1.128 and t = 

26.45 cycles, whereas, the envelope of the k = 10 curve falls below 
the exponential behavior due to numerical losses off the high- 


momenta end of the grid. 
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We now wish to examine steady- state distribution functions 
produced by the magnetic -pumping process. The first case we inves- 
tigate is one with the diffusion coefficient constant with momentum 
(Figure 6), as in our previous models. Particles are produced at a 
constant rate with low momenta values and are accelerated, on the 
average, toward higher momenta. Since the diffusion coefficient is 
constant, with the use of particle flux conservation, we may predict 
the steady- state shape of the distribution function. The, total 
momentum of a distribution of particles is described by Equa- 
tion ( 2 6 ), where, since the diffusion coefficient is not a function 
of momentum, f is momentum independent. Conservation of particles 
insures that, as an element of distribution between the low momentum 

values of P and P >- AP moves to higher momentum (and spreads in 
o o o 

momentum, into F to P + AP), 


f(P )AP = f (P)AP 
o o 


Using Equation (26) we see that the low-momentum boundaries of the 
element move to 


P -♦ P = P e 
o o 
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P + AP -* P + AP = 
o o 



+ AP o )e t//T 


(29) 


giving 

AP = AP e t//r . (30) 

o 

Inserting Equation (30) into the particle flux conservation relation 
(27) and using Equation (26) to eliminate the exponential, yields 

f(P )P 

f(P) = — |-a . (31) 

Thus, the distribution function should be of the form 

f(P) « i . (32) 

- , 903 

The distribution function of Figure 6, with a P falloff, is in 

excellent agreement' with this . 

The next cases considered concern diffusion coefficients &(P) 
that vary linearly with momentum. Figure 7(a) displays the results 
of a coefficient variation which produces less efficient energiza- 
tion at higher momenta than it does at lower momenta. As is evi- 
denced in the figure, we may expect particles produced at 
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low-momenta to "pile-up" as they are accelerated to higher momenta, 
causing the distribution to fall off slower than l/P. One inter- 
esting consequence of such a diffusion coefficient is the creation 
of distribution functions with positive slopes. In Figure 7(b) the 
diffusion coefficient varies in the opposite sense, with the energi' 
zation being more efficient at higher momenta than at lower. Here, 
as intuition would indicate, the distribution falls off faster than 
l/P. 

The final cases of momentum- dependent pitch-angle diffusion 
coefficients that we choose to examine concern coefficients that 
vary strongly (or discontinuously) at one momentum value, producing 
abrupt transitions between (momentum) regions of differing energiza' 
tion efficiencies. The steady-state distribution function of 
Figure 7(c) bridges the transition between a region of high energi- 
zation efficiency (at low momenta) and a region of considerably 
lower efficiency (at high momenta), with only a small effect at the 
transition momentum. If the diffusion coefficient is such that the 
regions of energization efficiency are reversed (Figure 7(d)), with 
low efficiency at low momenta and a sharp transition to high effi- 
ciency at high momenta, then a steady- state distribution function 
that may be interpreted as having "two temperatures" will result. 
Numerical simulations show that the difference in efficiency across 
the transition may be small yet still produce a noticeable break in 
the distribution at that momentum value. "Two- temperature" 
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distributions are not limited to the steady-state case, but will 
also evolve from an initial distribution function if particle pro- 
duction is not present. 

In this section we have compared (favorably) numerical simu- 
lations with theoretical predictions for the Alfvlnic model of 
magnetic pumping and have completed elementary exercises for the 
continuous-diffusion magnetic -pumping model, finding a maximum 
energization occurring at a "resonance" between the diffusion- and 
magnetic cycle -times, and explored the distribution functions 
generated by the process when the pitch- angle diffusion rate is 
allowed to be a function of particle momentum. We now apply the 
continuous -diffusion magnetic-pumping mechanism to particles trapped 
in the corotating magnetosphere of Jupiter. 
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3. PARTICLES IE THE JOVIAN MAGNETOSPHERE 

Due to the azimuthal asymmetry of the Jovian, magnetic field, 
particles -will drift around the planet through regions of differing 
magnetic field strength, and due to the presence of plasma waves, 
those particles will be pitch-angle scattered as they drift. Thus, 
for particle distributions trapped in the magnetosphere, there exist 
mechanisms both for changing the magnetic field strength (adiabatic 
motion) and for pitch-angle scattering (diffusion). 

We first determine the regions (if any) of the magnetosphere 
in which the energizing mechanism may be applicable. In order to 
model the equatorial magnetic field strength and to calculate the 
eorotational electrical potentials of Jupiter, we choose [as in 
Goertz, 1978] the Barish-Smith [1975] magnetic field model for the 
compressed dayside magnetosphere and the Goertz et al. [1976] model 
for the nightside magnetotail. (See Figure 8.) It is assumed that 
for distances less than 20 Rj (Rj « 7*10^ cm) the Jovian magnetic 
field is approximately dipolar, whereas the nightside model falls 
considerably below the proper dipolar values close to the planet. 
Hence, we shall apply our model only in the outer magnetosphere 
(beyond r = 20 R_) and assume the magnetotail model matches a dipole 
field somewhere within 20 Rj- Although these field models are 
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idealized we believe that they represent the important features of 
the asymmetric Jovian magnetosphere. 

The equatorial electric potentials (Figure 9) were obtained 
by integrating the corotational electric fields from r = » to the 
position of interest, that is 

V = J* ^ (v X S)dr , (33) 

CO 


where we will assume that the magnetosphere corotates in the region 
of interest (20 - 50 Rj). For the magnetotail model, integral (33) 
is written 


V = 
TAIL 



r 2 ' 7 


dr + constant 


(34) 


where the constant shifts the potential to match, at r = 20 Rj, the 
corotational electric potential of a dipole. Integrating, we find 
(r in cm) 


V = — 
TAIL c 




£ M 

C 20 Rj 


(35) 


where, 
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-k , 

a » 1.77 X 10 /sec 

5 30 3 

M b 4.225 gauss Rj = 1.449 x Hr gauss era 

M' b .65 gauss R 2 *^ = 2.4 x 10 2 ^ gauss era 2 ’ 1 . (56) 

The dayside corotational electric potentials are more difficult to 
obtain since the Barish-Smith magnetic field model equations must be 
solved by iteration at every position in space, hence the equatorial 
electric potentials were obtained by numerical integration (and 
again, the potential was matched to dipolar values at r = 20 Rj) . 

As charged particles drift in the magnetosphere, their total 
energy (kinetic plus electric potential) is conserved. For low 
energy particles (with much less kinetic energy than electrical 
potential), the drift is along equipotential surfaces. Hence, as a 
function of frontside radii, we may calculate the tailside drift 
radii and the front- and tailside equatorial magnetic field values. 
Knowing the front (maximum) and tail (minimum) magnetic field 
strengths gives us the ratio k (relation (25)), as a function of r 
(frontside), plotted in Figure 10. k.’th k(r), we are in a position 
to claculate the "energization time" of a distribution as a function 
of the frontside radius. We define the energization time t,-, to 'be 
the time, in number of Jovian rotations, that it takes to double the 
mean (nonrelativistic) energy of a distribution of trapped parti- 
cles. Numerically simulated values of t_ for the continuous 

Jli 
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diffusion model with various (ce -independent) diffusion coefficients 
are plotted as functions of dayside radius in Figure 11. Energiza- 
tion times for the Alfvenic model (giving the fastest possible 
energization for magnetic pumping) are plotted for comparison. 

To determine the regions of the magnetosphere in which the 
pumping mechanism energizes particles, we must compare the energiza- 
tion time t e to particle loss times. In the outer magnetosphere 
loss cone escape of particles is negligible [Coroniti, 1975]. We 
use the following value for the "normal" radial diffusion coeffi- 
cient in the outer magnetosphere [Coroniti, 1975] 

D I = 2 x 10" 10 1 ? / sec , (37) 


to estimate the radial diffusion loss time (see discussion below), 


T 


RD = 


1 SB 1 Bf 
B 3r LL f 3r 


1 BB 
B Br 


D — — 

LL 1 R 



(38) 


The result is the curve in Figure 11. Particles are energized in 
regions where the energization time t e is less than the radial 
diffusion loss time T™.. Note that the disclosure of energizing 
regions is very sensitive to the estimate of the radial diffusion 
coefficient, which, if estimated to be a factor of ten higher than 
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the value assumed, will lead to the prediction of no particle ener- 
gization in the 20 - 50 Rj region examined. 

Energy conservation for a drifting particle yields an upper 
limit to the kinetic energy that a particle may gain through the 
magnetic -pumping process. A low energy particle will f x I drift on 
an orbit which lies on an electrical equipotential surface, the 
drift carrying the particle through regions of differing magnetic 
field strength, hence the possibility of it being magnetically 
"pumped". A high energy particle will drift on a surface of con- 
stant magnetic field strength, eliminating the possibility of 
further magnetic pumping. To estimate this kinetic energy limit, we 
calculate the difference in electrical potential, AV, along a drift 
orbit which lies on a surface of a constant magnetic field and 
assume that the magnetic -pumping process will give energy to a 
particle as long as it has a kinetic energy 

W < q. AV . (59) 

For all B = constant drift orbits in the region from 25 Rj to 45 Rj, 
we find that AV — 2 MeV, so, for the magnetic field models chosen, 
the magnetic-pumping mechanism may, in that region, energize parti- 
cles (ions and electrons) up to energies of at least 1 MeV. 

An increase in kinetic energy due to the magnetic-pumping 
process will correspond to a radial transport of particles. The 
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total energy (kinetic plus electrical potential) being constant, the 
change in kinetic energy is related to a change in the potential 
energy by 

AW = - qE Ar . (4o) 

Using E = - Cl r B/c for the corotational electric field, and divid- 
ing by Ac, relation (40) is rewritten 

5W q il r B Sr /i,-\ 

at = c at > w 

giving the radial velocity of the transport, 


_ c aw 
V r - q r B at 


(42) 


2 

Knowing W = P /2m, we may write Equation (26) in terms of an energi- 
zation time, t e , as 


V t e 

W = W e 
o 


(43) 


enabling us to replace the derivative in Equation(42) , yielding the 
radial transport velocity 
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c W 

q ft r B T e 


(44) 


Note that, for energization, this corresponds to an outward (toward 
larger r) flow of ions and an inward flow of electrons. 

If we assume that, for any particle in a distribution func- 
tion, the chances of being energized or de-energize i during a cycle 
of the process are equal, then the radial displacements Ar of a par- 
ticle have zero mean and the radial transport may be described as a 
diffusion process. To approximate the radial diffusion coefficient 
we write 



ii&fl 

T« 


(45) 


where AL — Rj Ar and t r = 3*57 X 10 sec is the Jovian rotational 
period. We obtain AL from relation (40) and approximate the change 
in kinetic energy during one cycle by AW as 3W/dt t r with Equa- 
tion (43), yielding 



W 2 


2 2 2 2 
q E R t T, 


JB 


( 46 ) 


or, replacing the corotational electric field E, 
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_w \ 

r V 


2t r C 


2 2 2 
q 2 R J 


( 47 ) 


This type of radial diffusion preserves the total energy but not the 
first or second adiabatic invariants. 

Our treatment neglects the displacement of particles in the 
radial direction which is a necessary consequence of the pumping 
process. However, if the displacement AL in one rotation period is 
small (AL/L « l) our results should accurately describe the local 
change of characteristic energy. Particles can also diffuse 
radially by flux tube interchange motion or "normal" L shell diffu- 
sion conserving the first and second adiabatic invariant. If this 
process is faster than the L shell diffusion at constant total 
energy the adiabatic deceleration of outwardly diffusing particles 
will be greater than the energy gain by magnetic pumping. The char- 
acteristic energy of the plasma ions in the Jovian magnetosphere 
will decrease with increasing distance in those regions where L 
shell diffusion at constant first and second adiabatic invariant is 
faster than the L shell diffusion at constant total energy. "Nor- 
mal" L shell diffusion can be viewed, as a loss of particles at con- 
stant energy (for outward diffusion and a spectrum that decreases 
with energy), and its strength can be characterized, by a loss time 
given by Equation (38) . 
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For completeness, we look briefly at the pitch-angle distri- 
bution functions of particles (ions and electrons) corotating with 
the Jovian magnetosphere. The cases we choose to display concern 

k = B /B =4 orbits where, in one case the particles experi- 
°MAX °MIN 

ence continuous diffusion (Figure 12(a)), and in the other case 
experience increased amounts of diffusion in differing regions of 
their orbits (Figure 12(b)). All cases are associated with pitch- 
angle -independent diffusion coefficients. By comparison with 
measured Jovian pitch-angle distributions, these computer-generated 
distributions may be useful for estimating the absolute amounts 
(Figure 12(a)) of and spatial (local time) dependencies 
(Figure 12(b)) of the actual diffusion coefficients. 

In order to predict distribution functions in the Jovian 
magnetosphere, we construct a more realistic (pitch- angle and 
momentum- dependent) diffusion coefficient, describing whistler-mode 
scattering of electrons and ion cyclotron scattering of ions. The 
pitch-angle dependence is modeled after the terrestrial electron 
bounce -averaged diffusion coefficient of Lyons, Thorne, and Kennel 
[1972, Figure 4], where the diffusion of large pitch-angle particles 
is due to Landau resonance with the waves and diffusion of small 
pitch-angle particles is due to cyclotron harmonic resonances. To 
this diffusion coefficient we add a smooth low- momentum cutoff of 
the cyclotron resonance contribution (see Figure 2, top, of Lyons, 
Thorne, and Kennel) . To estimate the electron and ion velocities 
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at this cutoff, we use the Doppler- shifted cyclotron-resonance 
condition 



m 


where 

ua,k = frequency and wave number of wave 

ft = cyclotron frequency for relevant species 
v = particle velocity 



N = order of resonance 


Choosing the whistler-mode frequency to be » = l/lO ft and the ion- 
cyclotron-mode frequency to be in = l/lO ft^, we find, for N = - 1, 
the critical velocities to be, for electrons, 

l.l»fio B 

v c “ </kiT m„ ^/n~ * 

e o 

and for ions, 


( 49 ) 


(50) 


V 

c 


11 v ALFV^N = 11 Jkrrn Am ' 

v o p 


■where, 

n Q = plasma density 

m = mass of electron, proton 

e,P 

A = atomic number of dominant ion species 

In the vicinity of 50 Rj (B = 2 x 10 ^ gauss), we use the plasma 

-3 -3 

density limits of Barbosa et ,al. [19793 > = .035 cm - .14 cm , 

and find that the resonance cutoffs correspond to electron kinetic 
energies of from 1 KeV to 4 KeV and to ion kinetic energies of from 
8 KeV to 35 KeV, regardless of the species of the dominant ions. 

For both electrons and ions, the cutoff energies vary as B / n . Our 
diffusion coefficient is displayed in Figure 13. 

The distribution functions, both steady-state and time 
evolving, are of the "two-temperature" type with break points 
occurring near the cyclotron-resonance cutoff. A steady-state case 
is displayed in Figure l4(a) where particles are produced with 
momentum values below the cyclotron resonance cutoff transition 
region and are accelerated (on the average) to higher momenta. The 
distribution function maintains this "two-temperature" character 
through the extent of its orbital travel. This "two- temperature" 
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trait is omnidirectional; i.e., examination of the distribution 
function (f vs p) for any pitch angle will yield a similar shape. 
Figure l4(b) concerns the time development of a distribution under 
the action of our wave-particle diffusion coefficient. The double- 
temperature character sets in promptly (more rapidly if the diffu- 
sion is increased). In this case, as in the steady-state case, the 
distribution function’s traits are omnidirectional. 
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4. SUMMARY AND DISCUSSION 

We have extended the original model of Goertz [1978] to 
treat, under more realistic assumptions, the evolution of particle 
distribution functions in the Jovian magnetosphere. The model of 
Goertz assumed pitch-angle diffusion to occur only in certain 
regions of the Jovian magnetosphere (i.e., noon and midnight) and 
corresponds to the old model of Alfvln [1954] . In reality pitch- 
angle diffusion must take place in all regions of the magnetosphere, 
although its strength may be local-time dependent. We have numeri- 
cally followed the evolution of an arbitrary distribution function 
assuming that all particles stay on one drift shell. This is not 
exactly true but a reasonable approximation for energies below 
1 MeV. In the magnetic -pumping process the characteristic energy 
(mean energy) of the distribution increases exponentially with time 
and we have investigated the time constant t e of this increase as a 
function of magnetic field topology and pitch- angle diffusion 
strength. It should be stressed that as the mean energy increases 
the particles spread in real space, the radial diffusion being 
directly related to the pitch- angle diffusion strength and the field 
topology via t™* The mechanism does therefore not increase the 
total energy content of a flux tube if the plasma distribution is 
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uniform in space. However, if particles are injected in the inner 
region of Jupiter, the characteristic energy and energy distribution 
of the plasma varies with radial distance. If t e is smaller than 
the time for radial diffusion at constant first and second adiabatic 
invariants, the characteristic energy for ions increases with dis- 
tance, and for electrons, decreases with distance. Otherwise it 
decreases for ions and increases for electrons, but not as fast as 
expected on the basis of conservation of the first two adiabatic 
invariants. We will discuss this point in more detail in a future 
paper . 

We now will discuss the global aspects of our results. It is 
believed [Goertz, 1979; Hill, 1979] that Io is the main source of 
plasma in the Jovian magnetosphere, injecting ions at a considerable 
rate (10 ... 10 7 ions/sec). These maintain the "torus". How- 

ever, fluctuating fields will allow the torus ions to diffuse 
radially outwards (and inwards). Due to the azimuthal symmetry of 
the magnetosphere, within 20 Rj the magnetic-pumping time scale t e 
is very long and radial diffusion should preserve the first and. 
second adiabatic invariants. The plasma temperature should thus 
decrease with radial distance out to about 20 Rj, although not as 

_3 

fast as r . Beyond 20 Rj the magnetic pumping becomes effective 
d.ue to the azimuthal distortion of the magnetic field and the char- 
acteristic energy of the ions should increase with radial distance. 
The energy spectrum of charged particles should, be as shown in 
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Figure l4(a,b) with the break point occurring near 20 keV at 50 Rj. 
This spectrum is in agreement with that suggested by Belcher et al. 
[I98O]. We regard Figure l4(a,b) as the major new result of this 
work. 

A careful reader will have noticed that the magnetic pumping 
energizes only inward diffusion electrons. Magnetic pumping causes 
ions to diffuse outwards and electrons inwards. This constitutes a 
radial current which has to be closed if a steady- state rotation is 
to be maintained. We suggest that the current is closed by field- 
aligned currents carried by electrons out of the Jovian ionosphere. 
This is also implied in the model of Hill [1979] and has been sug- 
gested by us before [Goertz, 1979]* The ionospheric electrons can 
be trapped in the equatorial current sheet [Goertz, 1973, 1976] and 
provide a source for inwardly diffusing electrons which are strongly 
heated as they are transported. We will elaborate on these aspects 
in a future paper. In Figure 15 we give an overview of this model 
and sketch the expected distributions for electrons and ions at 
various locations in the magnetosphere. 
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APPENDIX I 


NUMERICAL TREATMENT OP ADIABATIC MOTION 
AND DIFFUSION 


Using the results of Section 1, we may follow the adiabatic 
motion through phase space of afty particle of a distribution trapped 
in a time-varying magnetic flux tube. An element of distribution at 

the old coordinates, (P°^,Q! 0LI) ) , is assigned to the new coordinates 
( p NEW^ Q NEW) ^ that is ^ n ^ p NEW^ Q NEW) _ n ( p OLD^ a OLDj ^ Numerically, we 

represent the old distribution function n(P 0L ^,oP L ^) on a finite 
grid in (p,0£ o ) phase space, assigning a certain fraction of the dis- 
tribution to each grid point. After the magnetic field is changed, 
we will find that, in general, unlike the old phase space coordi- 
nates, the new coordinates (pEEW^Q^NEW^ will not fall on grid points. 

NEW NEW\ 

Hence, a method of assigning the distribution n(P ,ot Q ) to the 
grid points near the position (pEEW^^NEW^ ^ obtained. The 

scheme used divides the distribution between the four grid points 
forming the rectangle containing the point (p^EW^^NEW^ . 
pNEW _ + v qjNEW _ + u ( see Figure l6 ) . The simple scheme 


chosen is 
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n(P c ,cg . (1 - g)(l - h^P^cf") 
n( p 0 ,Vi) “ e(l - h)n(P NEW ,o^ EW ) 
n(p c+l' a c> = (1 ' g)h ntP^aP*) 

n(P c+l^ a c+l ) = 8 h n(P NEW ,o^ EW ) (A.l) 

where g = g(u) and h = h(v). Note that this division conserves the 
number of particles 


n(P c ,0!c) + n(P ft ,a 


c * c+1 


+ n < P = + l> a o + l) 


) + -(r ctl ,“o) 

= n(F« of™) 


(A.2) 


We will choose g(u) and h(v) such that scheme (A.l) conserves two 

2 

further quantitites, P and Q^. These choices were found to give 

the greatest accuracy in reproducing a transformed-untransformed 

2 

distribution function. Solving the P conservation equation 


pNEW 2 ^pNEW ojNEW^ _ p 2 ( p a j n (p a ) 

s ’ o ' c ' c ' cr c c’ C4-1' 


.2 


P c+1 n< ' P c+l ,a c^ + P c+1 n ( p c+l ,a c+l^ 


(A.5) 


and the a conservation equation 


4o 


<P n( P “ -C n(p,a) n(p ,a ,) 

o ’o' c ' c* c' c+1 ' c # c+1' 

+ “c n < P c + l’ a c> + “o + l n < P c + l> a = + l> 
with the insertion of relations (A.l) yield, respectively, 


h = 


pNEW 2 _ p 2 

_c 

2 2 

P - P 
c+1 c 


(A.5) 


and 


g = 


c^-a 
o c 

a . - a 
c+1 c 


(A.6) 


or 


g = u 


(A.7) 


Another method may be used to follow the adiabatic transfor- 
mation in phase space of the distribution function. At each grid 
point of the new coordinates (P ,Qf , , we look back m time on 


the old (? 


OLD OLD 


) grid to ask where the element of distribution 


^pNEW^QjNEWj carne f romt Unlike (p NEW ,oP 3W ), the source point 


4i 


(pOLD^OLD^ not lie on a grid point, and n(P°^,a^' D ) must be 
estimated by an interpolation. Interpolation schemes using up to 
34 points on the (P 0LD ,a^ L ^) grid were tried, and two drawbacks were 
found to make the method unacceptable. The distribution function 
was not conserved in the transformation and, more hindering, the 
magnetic field had to be changed slowly in time in order to prevent 
large shape changes in the distribution function which introduced 
errors in the interpolation scheme. 

For numerical treatment of the pitch-angle diffus~'?n Equa- 
tion (22), we center it in time (implicit scheme), writing a set of 
equations 


N+l 


n. 


N 


- n. 


At 


da 


o <- 


dn 

da 


■jW+l/2 


o J 


da 


[H ni: 


N+l/2 


(A.8) 


where the superscripts (subscripts) represent the time (pitch -angle) 
grid points. The pitch-angle derivatives will be represented as 


da 


\ SsJ 

_i. r 

30 da 

o J 

Aa [_ 


1+1/2 da 


5n i -l/2 ' 


i-l /2 da 


^ Cn i + i ( *i + i + V + V- *i + i 


" 2j8 i " + + ^i-i^ > ( A -9) 
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and 


A 3H, an, 

doi~ n ^i = IcT n i + H i W 
O 0 o 


5*5 tn i + i H i + V H i + i - H i-i> 

+ H i» • 


(A. 10) 


Centering these in time about the point N + l/2 (we use n ^ +1 /2 = 
l/2 n^ + ^ + l/2 ri^), and inserting them into Equation (A. 8) yields 


O- K ^i + i “i)-*® + “H 1 + + 2i i + 


JI+lr 


- H A«(H. . - H. )} + n^{- k(j 8 . + fi. .) + K ha H.} 
' l+l i-l /J l-l 1 - ' l 1-1' i J 


■ n Ll tK(4l i + l + V + “ H i 3 + "i' 1 + “ ( - *1 + 1 


- 2 1 0 i - jD i _ 1 } + + fi.^) - H Aa H.} , (A. 11) 


where we have defined 


h = 


At 
4(Aa)‘ 


(A. 12) 


The implicit scheme (A.ll) may be written in tridiagonal form 
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N+l , . N+l N+l 

a. n. . + b. n. + c. n. _ = d. 

1 i+l li l i-I i 


(A.13) 


where 


a = C. - k Aa H. 
1 i 


b = 1 + C 2 - » ZW(H 1+1 - H ± _ 1 ) 


e = C, + k 40 H. 
3 i 


W w ro- 
ds t n n. + t 0 n. + t, n. , , 

1 i+l 2 l 3 i-I 9 


(A.l4) 


and 


c i = • K(<, i+i + V 

°£ , *(« 1+1 4 2 i i + 

C = - »(*! + 


t, = 


-CL + k AO H. 

1 l 


tg = 1 - C 2 4 K 4a(H 1+ , - H^_ x ) 


t, = 


- C, - k Aa H. 
3 i 


(A.15) 


Before looking at the algorithm used to solve the diffusion 


equation, we examine the necessary boundary conditions. The 
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requirement that the total number of particles in the flux tube be 
constant may be expressed 


tt/2 Bn(a ) 

J* — da 

j at o 

o 


= 0 


(A.l6) 


or, using Equation (22), as 


a 


Bn 

aa 


,tt/2 


+ [H n] 


a =o 

o 


tt/2 

a =o 

o 


(A. 17) 


This we may satisfy by requiring a(a ) and n(a Q ) to vanish at the 
boundaries, a^ = 0, tt/2. Note that n(a Q ) being zero for pitch 
angles 0 and tt/2 is consistent with the expression (l8) relating 
n(a ) to the distribution function f (a ) . To insure that n(a ) 
remain zero on the boundaries for all times, we constrain Bn(o)/Bt 
and Bn(n/2)/at to be zero. Applying Equation (22) at the two 
boundaries, we see that 


Bn 

aa aa 

o 



BH 

aa 

o 


n + H 


3n 

aa 


a 

o 



) 


. (A.18) 
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This is satisfied by insisting that d£/dQ! o , ^( a 0 ) (already 
required), n ( a 0 ) (already required), and H(c< o ) (consistent with 
definition ( 21 )) all vanish at the boundaries ct Q = 0 and a Q = tt/2. 

We now examine the simple tridiagonal algorithm used to solve 
the set of diffusion equations (A. 15 ) (see, for example, Roache 
[1972], Appendix A), consistent with the above boundary conditions. 
We assume there exist solutions of the form 


N+l 

n. 

x 


e. 

1 


N+l 

i +1 


+ f. 
x 


(A. 19 ) 


where we note that, here, f. is not the distribution function. We 
use this assumption at the (i - l)th spatial point to replace n^*^ 
in the tridiagonal equation (A. 15), resulting in 


N+l /. \ N+l „ 

n. , , + (b. + c. e. . )n. + c. f . . 

X+l v X X X-l' X X X-l X 


= d. 


(A.20) 


Comparing this to the assumption (A. 19 ), we find the relations 


e. 

x 



a. 

x 


+ c. e. . 
x x-l 


and 


f. 

x 



b. + c. 
x x 



(A. 21) 


If e n and f. (cc =0 boundary) are known, all other e. and f. values 
1 1 o ' xx 

(Recall that all a., b., c., and d. values are 

' X* X* x’ X 


may be generated. 


known. ) The boundary condition on n(a Q ) at ot Q = 0 is n(0) = 0, 
which, with definition (1.19), requires 

e l = ° , = 0 . (A. 22) 

Using these values for starting points, we calculate all e i and f^ 
values for i=2toi=M-l, where M is the number of grid points 
used to represent the pitch- angle coordinate (i = l<*0! =0; 
i = M * ot Q = tt/ 2). The boundary condition on n(a o ) at a Q = tt/ 2 is 
n(Tr/2) = 0, written as 

4 +1 = ° • (A.23) 

Using relation (A.19) we may now generate, since all f^ and e^ 

venues are known, all n^ values, from i = M - 1 to i = 2. Thus, for 

W+l 

time step N + 1, all n^ values are obtainable. 

We test our numerical pitch- angle diffusion scheme by follow- 
ing the evolution in time of an initially anisotropic distribution 
function (df/ciC^ ^ 0) . Figure 17 shows that after 100 time steps 
the distribution function is indistinguishable from isotropy. 
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APPENDIX II 


LIMIT OP CONSTANT ISOTROPY 


Our aim is to show that, for the model similar to that of 
Murty and Varma, in the limit of forcing the distribution to remain 
isotropic (infinite continuous diffusion) throughout the magnetic 
field cycle, there is no net momentum gain. To motivate this, we 
first examine the case of isotropizing at four stages (as opposed to 
two stages in the Alfvenic model) per cycle. The magnetic field 
values (equatorial) of the cycle will be 


B 


B, 


B„ 


V B 1 


B 


(A.24) 


Conserving the first adiabatic invariant as the magnetic field 
changes from B q to B^, we find 



1 

3 



o 


(A.25) 


and upon isotropizing, find 
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(A. 26) 


The adiabatic change in field strength from to Bg yields 


p l - i ( 2 r + 1 )( 2 r + x ) - 

3 o 1 


(A. 27) 


and upon isotropizing, 


P = 


P“ = 


«2 


I J p o ( 2 r - p )( 2 r + x ) 
H p c( 2 B , + 1 X 2 r + 1 


(A. 28) 


The field change from Bg to B^ gives 


B, 


B, 


p 3 ■ 3 p o ( 2 T * P X 2 T + p )( 2 ^ + 1 


(A. 29) 


Similarly, isotropizing and the field change from B^ to B^ will 
yield 
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-= \ p 2 ( 2 F + ^ ^ + "X 2 X + X X 2 X * ^ ' <A ' 30> 


Assuming that B $ = B ] _ and B 4 = B q , we 
tion per cycle to he 


find the fractional moment iza- 



(A.31) 


If the distribution is isotropized at N steps during the 
cycle, rather than at four, Equation (A.3l) becomes 



N/2 

tt ' r 


B. 

i 


l 

1 


5 + 2 (.hT: + 

i-l 


\=1' 

V- 


1/2 


(A.32) 


where we have assumed B^ 2 j + 1 = B (N/2)-i* For sim P llc:Lty > we let 
the magnetic field strength change linearly with time from B = B q to 

B = \/2’ 811(1 then decrease l ineaxly back from B = B N/2 to 

■r - B = B . Thus, we may write (for i £ N /2) 

D N o 


B. = B + i AB , 

i o 


(A.33) 


where 
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AB = 


B M/2 ~ B o 

N 

2 


with this, we may substitute 


B. B. . B. B. - AB , atA 2 

i_ , i-l l , _i a. p , ( ABN 

b7“ + B. = B. - AB + B. " 2 + \B. 

l-l ii l l 


in Equation (A. 32) to get 

^»/ 2 . 


l 

P o 3 S/£ 


1 

1 


5 + 4+2 

\iJ . / 
1 


2-i 


, 1/2 


Replacing B i by Bj^. = B q (which will give an upper 
hence to P N / P 0 )> v e may write 


■TSS ^ 1 


P o 3 H / 2 


_ _ ( AB\ 2 

n/4 

2 /ABN 2 ' 

9 + 2 (b ) 

*- O 

Si 

1 + 9 'sB ) 

L o J 


,N/1 


Using the binomial expansion, we rewrite this as 



(A. 34) 


(A. 35) 


(A. 36) 


bound to AB/B^, 


(A.37) 


(A. 38) 
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Replacing AB with its definition, relation (A. 34), the fractional 
momentum gain per cycle is 

i + (i + 2k + k 2 ) , (A. 39) 

o y 

where k = b jv[ax/ N = ^/2^o ( E 9 uation ( 2 3) ) • In the limit of 
N -» « (forcing constant isotropy), we find 



P 

o 


-+ 1 


(a.4o) 


resulting in no net moment™ gain over the cycle. 
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FIGURE CAPTIONS 

2 

Figure 1. The magnetic field model, B = B q + l/2 a z . Shown for 

B /B =10. Note that the two fields are equal at 

°MAX °MIN 

z = z . 
c 

Figure 2. Comparison of computer-simulated results with analytical 
values for the Alfvenic model. P = momentum after one cycle, 

P = initial momentum, k = B / B . Simulation results 

° °MAX °MIN 

for the continuous-diffusion model, with the diffusion coef- 
ficient chosen to give a maximum energy gain, are also 
displayed. 

Figure 3. The number of particles in a flux tube, n(P,a ), in the 
Alfvlnic magnetic-pumping model, (a) Initial isotropic dis- 
tribution, (b) distribution after adiabatic motion due to 
increase in magnetic field strength, (c) distribution is 

pitch- angle diffused to isotropy, holding B at B , 

° °MAX 

(d) distribution after adiabatic motion due to decrease in 
magnetic field. 


Figure 4. Computational results of continuous diffusion model. 

P = momentum after one cycle, P Q = initial momentum. For the 
parameters of our numerical treatment £ = 0.282 gives 

t dif = t cycle* 

Figure 5- Time dependence of mean momentum of distribution in the 
continuous diffusion model. 

Figure 6. Steady- state distribution function for a momentum- 
independent diffusion coefficient. 

Figure 7- Steady-state distribution functions for momentum- 

dependent diffusion coefficients, (a) Coefficient linearly 
decreases such that the energization process is less effi- 
cient at higher momenta, (b) Coefficient linearly increases 
such that the energization process is more efficient at 
higher momenta, (c) Diffusion coefficient with discontinuity 
such that the energization process is more efficient at lower 
momenta. Note irregularity in vicinity of discontinuity. 

(d) Diffusion coefficient with discontinuity such that the 
energization process is more efficient at higher momenta. 
Definite two- temperature appearance. (Note (d) is semi-log.) 

Figure 8. Models used to approximate the Jovian magnetosphere. 
Barish-Smith for dayside, Goertz et al. for magnetotail. 
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Figure 9 . Corotational electrical potentials in the Jovian mag- 
netosphere. Calculated using the magnetic field model of 
Figure 8 . 

Figure 10. Magnetic ratio k = B /B in the Jovian magneto- 

°MAX °MIN 

sphere as a function of the dayside radius. 

Figure 11. Energization times T„ (for various amounts of diffusion) 

ill 

and radial diffusion loss times t rd as functions of the 
Jovian dayside radius. Energization is predicted in regions 
■where t e is less than t r:d . 

Figure 12. Jovian pitch-angle distributions, (a) for various 

amounts of diffusion per orbit and (b) for increased diffu- 
sion in various magnetosphere regions. Pitch-angle inde- 
pendent diffusion coefficients were used. 

k = B /B =4. 

°MAX °MIN 

Figure 13. Pitch-angle diffusion coefficient modeled for electron 
resonances with whistler modes or ion resonances with ion- 


cyclotron modes. 
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Figure l4(a). Steady-state distribution functions for wave -particle 
diffusion coefficient. In the case of electrons interacting 
with whistler modes, the transition region corresponds to 
electron kinetic energies of from 1.1 to 5.5 KeV. In the 
case of ions interacting with ion-cyclotron modes, the 
transition region corresponds to ion (any species) kinetic 
energies of from 10.9 to 54.3 KeV. (b) Time development of 
distribution function using wave resonance diffusion 
coefficient . 

Figure 15. Schematic overview of outer magnetosphere with magnetic- 
pumping process transporting ions outward and electrons 
inward. 

Figure 1 6 . The momentum-pitch- angle phase-space grid. Distribution 
at is divided between the four solid points. 

Figure 17. Time action of the pitch-angle diffusion equation. Note 
the conformity of the large time result to the exact solu- 
tion, n(a) = (constant) sin a cos a T(a) corresponding to 
isotropy. 
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